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Abstract. Theoretical foundations of the problem of quantum spin tunnel- 
ing in magnetic nanostructures are presented. Several model problems are 
considered in detail, including recent new results on tunneling in antiferro- 
magnetic nanoparticles and topologically nontrivial magnetic structures in 
systems with reduced dimension. 



1. Introduction 

It is well known that magnetic ordering is an essentially quantum phe- 
nomenon. According to the Bohr - van Leeven theorem (see, e.g., [1]), 
the magnetization of a thermodynamically equilibrium classical system of 
charged particles is zero even in presence of an external magnetic field. Clas- 
sical theories of magnetic properties were based on certain assumptions go- 
ing beyond the limits of classical physics (e.g., the existence of stable micro- 
particles with nonzero magnetic moment assumed in Langevin's theory of 
paramagnetism [1]). The nature of magnetic ordering was revealed only 
after the discovery of modern quantum mechanics in the works of Heisen- 
berg, Frenkel and Dorfman. In 30s, many remarkable results were obtained 
within the microscopic quantum theory: Bloch [2] predicted the existence 
of magnons and low-temperature behavior of magnetization; Bethe [3] was 
able to construct the complete set of excited states for a spin-^ chain, in- 
cluding nonlinear soliton-type excitations (spin complexes). 

The 'undivided rule' of the quantum theory of magnetism lasted only 
till 1935, when in the well-known work Landau and Lifshitz [4] formulated 
the equation describing the dynamics of macroscopic magnetization of a 
ferromagnet (FM). When deriving the Landau-Lifshitz (LL) equation, a 



quantum picture of magnetic ordering was used, particularly, the exchange 
nature of spin interaction, but the LL equation itself has the form of a 
classical equation for the magnetization M. Later on the basis of the LL 
equation the macroscopic theory of magnetism was developed and enormous 
number of various phenomena were described [5, 6] (an overview of modern 
phenomenological theory of magnetically ordered media can be also found 
in this book in the lecture by V.G.Bar'yakhtar). 

This lecture presents an introduction to the foundations of a new, fast- 
developing topic in the physics of magnetism, Macroscopic Quantum Tun- 
nelling (MQT). Let us first address briefly the scope of problems belonging 
to this field. MQT problems can be roughly divided into two main types. 
First of all, there are phenomena connected with the underbarrier transition 
from a metastable state, corresponding to a local minimum of the magnet 
energy, to a stable one. Such effects were observed in low-temperature re- 
magnetization processes in small FM particles as well as in macroscopic 
samples (due to the tunneling depinning of domain walls), see the recent 
review [7]. Such phenomena of "quantum escape" are typical not only for 
magnets, e.g., quantum depinning of vortices contributes significantly to 
the energy losses in HTSC materials [8]. 

Here we will concentrate on another type of phenomena, the so-called 
coherent MQT. To illustrate their main feature, let us consider a small FM 
particle with the easy axis along the Oz direction. If the particle size is 
small enough (much less than the domain wall thickness Ao), the parti- 
cle is in a single-domain state, because the exchange interaction makes the 
appearance of a state with magnetic inhomogeneities energetically unfavor- 
able. Then, from the point of view of classical physics, the ground state of 
the particle is twofold degenerate. Those two states correspond to two local 
minima of the anisotropy energy and are macroscopically different since 
they have different values of macroscopic magnetization M = ±Moe z . The 
situation is the same as in the elementary mechanical problem of a particle 
in two-well potential U(x) having equivalent minima at x = ±a, see Fig. 
1. In classical mechanics the minimum of energy corresponds to a particle 
located in one of the two local minima of the potential. 

However, from quantum mechanics textbooks it is well known that the 
actual situation is qualitatively different: the particle is "spread" over two 
wells, and the ground state is nondegenerate [9]. One can expect that the 
same should be true for a FM particle: its correct ground state will be 
a superposition of "up" and "down" states, and the mean value of mag- 
netization will be zero. Such picture was first proposed by Chudnovsky 
[10]; further calculations showed [11] that such effects are possible for FM 
particles with rather large number of spins (about 10 3 10 4 ). The tunnel- 
ing effects, according to the theoretical estimates [12, 13], should be even 



more important for small particles of antiferromagnet (AFM); the effects 
of quantum coherence in AFM particles were observed in Ref. [14]. 

Thus, an important feature of quantum mechanics, a possibility of un- 
derbarrier transitions, can manifest itself in magnetic particles on a macro- 
scopic (strictly speaking, mesoscopic) scale. Maybe even more interesting 
is the manifestation of another characteristic feature of quantum physics, 
viz. the effects of quantum interference. Such effects arise in the problem of 
MQT in magnetic nanostructures and can partially or completely suppress 
tunneling, restoring the initial degeneracy of the ground state [15, 16]. We 
wish to remark that understanding that motion of particles along very dif- 
ferent classical trajectories can "sum up" in some sense and yield an inter- 
ference picture was one of the crucial points in the development of quantum 
mechanics, and a considerable part of the well-known discussion between 
Bohr and Einstein was devoted to this problem. Besides the importance of 
the tunneling phenomena in magnets from the fundamental point of view, 
they are potentially important for the future magnetic devices working on 
a nanoscale. 

In the present lecture we restrict ourselves to discussing the problems 
of coherent MQT in various mesoscopic magnetic structures. The paper is 
organized as follows: Sect. 2 contains the elementary description of the in- 
stanton formalism, traditionally used in the theoretical treatment of MQT 
problems. Since the instanton approach, though being the most straight- 
forward one, is based on rather complicated mathematical formalism, we 
will discuss it in parallel with simple and widely known semiclassical ap- 
proximation of quantum mechanics. The point is that those two approaches 
are equally adequate for treating the problem of MQT in small particles, 
and the "standard" semiclassical calculations, easily reproducible by any- 
body who learned foundations of quantum mechanics, may be helpful for 
understanding the structure of the results derived within the instanton 
technique. Further, in Sections 3 and 4 we discuss the problem of MQT 
in ferro- and antiferromagnetic small particles, with a special attention 
to the interference effects. For the description of AFM we use simple but 
adequate approach based on the equations for the dynamics of the antifer- 
romagnetism vector I. This approach easily allows one to keep trace of the 
actual magnetic symmetry of the crystal; the symmetry is lowered when ex- 
ternal magnetic field is applied or when certain weak interactions, e.g., the 
so-called Dzyaloshinskii-Moriya (DM) interaction, are taken into account, 
which leads to quite nontrivial interference phenomena. Section 5 is de- 
voted to the analysis of coherent MQT in "topological nanostructures," i.e. 
static inhomogeneous states of magnets with topologically nontrivial distri- 
bution of magnetization; among the examples considered there are domain 
walls in one-dimensional (ID) magnets [17, 18, 19], magnetic vortices [20] 



and disclinations [21] in 2D antiferromagnets, and antiferromagnetic rings 
with odd number of spins [22]. For those problems, when the description 
of tunneling involves multidimensional (space-time) instantons, there is no 
alternative to the instanton approach and its use is decisive. Finally, Sec- 
tion 6 contains a brief summary and discussion of several problems which 
are either left out of our consideration or unsolved. 

2. Basics of Tunneling: With and Without Instantons 

For the sake of the presentation completeness, let us recall briefly the main 
concepts of the instanton technique, since we will extensively use them 
below. 

In quantum field theory, the propagator, i.e., the amplitude of proba- 
bility P\2 of the transition from any given state with the field configuration 
^>a{x) at t = to another state <Pb{x) at t = to is determined by the path 
integral 



is the action functional. Here C is the Lagrangian density, and the inte- 
gration in (1) goes over all space-time field configurations (p(x,t) satisfying 
the boundary conditions (/?(x,0) = ^>a{x) and (p(x,to) = </?#(x). (We leave 
out the problem of a consistent definition of the measure T>(p that arises 
for systems with infinitely many degrees of freedom, keeping in mind that 
we are going to talk about the application of field theory to the physics of 
spin systems on a discrete lattice, and thus all necessary regularizations are 
provided by the lattice in a natural way.) 

Instead of working with the propagator (1) in usual Minkovsky's space- 
time, it is convenient to make the Wick rotation t — ► it (essentially this 
procedure is an analytical continuation in t), passing to the Euclidean space- 
time. Then one has the Euclidean propagator 



The main contribution to the path integral comes from the global minimum 
of the Euclidean action functional A enc \. This minimum corresponds to a 
trivial solution ip = ipo = const, where ip$ determines the minimal energy 
of the system. However, if several different values of 990 are possible, it is 
often important to take into account the contribution from the local minima 
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of the Euclidean action as well. Such a local minimum can correspond, 
e.g., to a trajectory ip = ^mstfV) connecting two possible 990 values; it is 
clear that the probability Pab wm contain the factor exp{— -4 e uci[Vinst]}/^- 
Such a contribution can be calculated in a semiclassical approximation and 
describes effects which cannot be accessed by means of the perturbation 
theory. 

We will illustrate the above arguments on the example of a simple 
quantum-mechanical problem. Consider the motion of a particle of mass 
m in a symmetric two-well potential U(x) of the type shown in Fig. 1, with 
two equivalent minima at x = ±a. Following the popular choice [23], we 
will assume this potential in the form 

U{x) = \{x 2 - a 2 ) 2 , (2) 

where the parameters A and a determine the height and width of the barrier 
between two wells. This model is described by the Lagrangian 

After passing to the imaginary time, the Euclidean action is easily obtained 
in the form 

And = £ dr{\m + U(x)} . (4) 

The classical (global) minimum of this functional is reached at x = a or 
x = —a. Equations of motion for the action (4) 

d 2 x _ dU 
dr 2 dx 

correspond to the particle moving in the potential —U(x), so that x = ±a 
are maxima of this effective potential, and there exist classical low-energy 
trajectories connecting them. Such trajectories represent local minima of 
the Euclidean action functional and are called instantons. They can be 
easily found in implicit form, 

i dx \mx)) =T - ro ' (5) 

where tq is an arbitrary parameter determining the "centre" of instanton 
solution. For many potentials the integration can be performed explicitly, 
e.g., in case of (2) one obtains 



x = ±atanh[a;o(T — To)/ 2] , 



(6) 




Figure 1. A two- well potential U(x) with equivalent wells at x = ±a. Semiclassical 
treatment of tunneling is possible if the amplitude of zero-point oscillations S <C a. 

where ujq = (8Aa 3 /m) is the frequency of linear oscillations around one 
of classical minima. Euclidean action for the instanton trajectory can be 
written as 

A = I j2mU(x) dx. (7) 

J —a 

For the model (2) one has Ao = 8a 3 V Am/3. Thus, instantons are very much 
like solitons with the difference that they are localized in time. Trajectories 
(6) begin at r — ► — oo in one of the minima of U (ip) and end at r — > +oo in 
the other one; the contribution of those trajectories is responsible for the 
tunneling splitting of the lowest energy level in the two-well potential. In- 
deed, the tunneling level splitting is proportional to the matrix element t\2 
of the transition from one well to the other, and the probability amplitude 
of such a transition is given by the path integral from x = a to x = —a. 
It is thus clear that the contribution of a single instanton to the transition 
amplitude is proportional to e~ A °/ h . 

The full calculation of this amplitude, however, is more complicated 
and should take into account not only the instanton trajectories but all 
trajectories close to them. Further, the full variety of multiinstanton paths 



which bring the particle from one well to the other should be taken into ac- 
count. If the problem is semiclassical, i.e. Ao/h is large and the probability 
of tunneling is small, integration over "close" trajectories can be described 
as an effect of small fluctuations above the instanton solution. Even this, 
usually elementary, problem of integrating over small (linear) fluctuations 
is nontrivial in case of instantons, because some of those fluctuations do 
not change the action. Particularly, from (6) it is easy to see that changing 
the position of instanton centre tq has no effect on Aq. Such "zero modes" 
always arise in instanton problems and their contribution requires a special 
analysis. Detailed description of this technique would take us out of the 
space limits, and we refer the interested reader to textbooks and review 
articles (see, e.g., [24, 23]). 

We will attempt to get the correct result for the probability amplitude 
Pab by means of the "traditional" quantum mechanics (without use of path 
integrals and instantons). First, let us note that, due to the symmetry of 
the potential U(—x) = U(x), two lowest levels correspond to even and odd 
eigenfunctions ip s (x) and ip a (x), with the energies E s and E a , respectively. 
Multiplying the Schrodinger equation for ifj s by ip a and vice versa, then 
taking the difference of those two equations and finally integrating over x 
from to 00, one obtains the relation 

= / Ipal/Js dx , (8) 
x=0 JO 

which is exact and is nothing but a mere consequence of the symmetry 
properties. 

It is natural to try to use a semiclassical approximation. The semiclas- 
sical result is given, e.g., in a popular textbook by Landau and Lifshitz 
([9], see the problem 3 after §50). According to that result, E s — E a = 
{Tiujq/tx) exp{— A' /fi}, where ujq = [k/m] 1 / 2 , k = (d 2 U/dx 2 ) x=a and A' = 
J^q, [2m(U(x) — E)] 1 ^ 2 , here a' is the turnover point of the classical tra- 
jectory with energy E (corresponding to a non-split level) defined by the 
equation U(a') = E. However, this result is not adequate for our problem, 
and it does not coincide with the result of instanton calculation. The point 
is that, surprisingly, the problem of tunneling from one classical ground 
state to another is not semiclassical: semiclassical approximation cannot 
be directly applied to the ground state wavefunction inside one well. 

Therefore we will do as follows: let us represent the wavefunctions inside 
the barrier region as symmetric and antisymmetric combinations of the 
WKB exponents, 



(E a — E s 



d^pa 



dx 



where \p\ = y / 2m[U(x) — E}. Those wavefunctions can be used inside the 
entire barrier region, except narrow intervals \x ± a\ < 5 near the well 
minima, where 5 = (h/muio) 1 / 2 is the amplitude of zero-point fluctuations. 

On the other hand, if the condition a 5 is satisfied, then for the 
description of the wavefunction inside the well any reasonable potential 
U(x) can be replaced by the parabolic one, U(x) — > (k/2)(x ± a) 2 . Then 
in "non-semiclassical" regions one may use well-known expression for the 
ground state wavefunction of a harmonic oscillator, 

(ir6 2 )- 1 / A exp[(x±a) 2 /25 2 ]. (10) 

Thus, in the regions a 2 S> (x ± a) 2 » 5 2 both the expressions (9) and 
(10) are valid. Then, normalization factors C s>a can be determined from 
the condition of matching (9) and (10) in the two above-mentioned regions, 
and after that the integration in (8) can be performed explicitly. After some 
amount of algebra the tunneling level splitting can be represented in the 
form 

/2 { f a ~ s I U"(a) ] f 1 1 

where the quantity Aq = J*® dx y / 2mU(x) coincides with the Euclidean 
action for the instanton trajectory. 

One can see that the difference between the formula (11) and the usual 
semiclassical result consists in the pre-exponential factor containing the 
integral of the type / dxU~ l l 2 {x). It is clear that the main contribution 
into this prefactor comes from the region x ~ a, where the integral can 
be approximated as f a ~ 6 dx/\a — x\, so that it diverges logarithmically at 
<5 — > 0. Thus for any potential U (x) the prefactor can be represented in the 
form C(a/§) or, equivalently, {CAq/K) 1 ^ 2 . Here C is a numerical constant 
of the order of unity, it can be easily calculated for any given potential 
U(x). So, finally we arrive at the following universal formula: 

For the model potentials U = X(x 2 — a 2 ) and U = 2Uq sin 2 x the value of 
C is equal to y/3 and y/2, respectively. 

The formulas (11,12) give the desired result for any two-well potential 
with sufficiently large barrier. The main feature of this result is the pres- 
ence of an exponentially small factor. The small parameter of the MQT 



problem is h/Ao, which can be represented as a ratio of the zero-point fluc- 
tuations amplitude to the distance between wells, (h/Ao) ~ (5/a) 2 . The 
expression e~ a ' is non-analytical in the small parameter, and thus the 
MQT phenomenon cannot be obtained in any order of the perturbation 
theory. We wish to emphasize that the correct result is roughly (Ao/h) 1 / 2 
times greater than that following from "naive" semiclassical formula. This 
large additional factor appears due to the contribution from the regions 
close to the minima of the potential, where the motion is not semiclassic. 
Let us try to understand this in the instanton language. 

As we mentioned before, the small exponential factor exp(— Ao/h) arises 
immediately in the instanton approach; the main problem is to compute the 
pre-exponential factor, which is determined by the integration over all small 
deviations from the instanton solution. Those deviations are of two types: 
real fluctuations of the instanton structure, which increase the Euclidean 
action, and "zero modes" which correspond to moving the instanton centre. 
It is rather clear that "nonzero" modes have a characteristic energy of the 
order of huo, and that the quantity ujq has nothing to do with the zero 
mode. Thus, it is obvious that the factor Tlujq arises from the integration 
over all "nonzero" modes, and the large factor (Ao/h) 1 / 2 arises due to the 
zero (in our case - translational) mode. Such a "separation" naturally arises 
in rigorous calculations [23, 24]. 

It is remarkable that the above result can be generalized to the case of 
much more complicated problems involving space-time instantons (which, 
as we will see later, is important for the problem of MQT in topological 
nanostructures). For any instanton all nonzero modes yield a factor like 
huJo, and each of the zero modes yields the factor (Ao/h) 1 / 2 [23, 24], so 
that the final result can be reconstructed practically without calculations 
(up to a numerical factor of the order of unity) . 

To illustrate one more feature typical for tunneling problems, let us 
consider another model [24]: a particle of mass m which can move along 
the circle of radius R, so that its coordinate is determined by a single 
angular variable <p, < (p < 2ir, in the two-well potential 

U(ip) = U (1- cos 2ip). (13) 

The model is described by the following Lagrangian: 

H rf (t) 2 -^>- < i4 > 

The classical Lagrangian can be modified by adding the arbitrary full 
derivative term, e.g., 

L»L + y%, (15) 



which of course does not change the corresponding classical equations of 
motion. However, adding the full derivative (15) changes the definition of 
the canonical momentum conjugate to ip, which, as one can easily check, 
leads to a considerable change in the Hamiltonian of the corresponding 
quantum-mechanical system after canonical quantization: for nonzero 7 the 
correct Hamiltonian would be 

"=^KW 2+C%) - (16) 

Thus, there is no one-to-one correspondence between classical and quantum- 
mechanical systems: several quantum systems can have the same classical 
system as a classical limit. 

For this model problem the instanton trajectories can be written down 
explicitly: 

cos if = <Ji tanh[cj(r — Tj)] , (17) 
u = (4[/ /mi? 2 ) 1/2 . 

where T{ is the arbitrary parameter determining the instanton position in 
the imaginary time axis and Cj = ±1 is the topological charge distinguishing 
instantons and antiinstantons; the instanton action is finite and is given by 
A = (8mR 2 U ) 1/2 . 

The importance of the full derivative term (15) can be most easily un- 
derstood in terms of instantons. Indeed, let us consider the tunneling ampli- 
tude P\2 from the <p = well to <p = tt one: it is clear that the contribution 
to this amplitude is made equally by instantons (with ip changing from to 
tt) and antiinstantons (with tp changing from to —tt). However, the term 
(15) becomes an imaginary part of the Euclidean action and leads to the 
additional factor e' ml / h associated with the instanton contribution and a 
similar factor e"" 7 / fi for antiinstanton paths. Thus, the resulting transition 
amplitude for nonzero 7 is modified as follows: 

P12 = [Pl2] 7 =0COs(7T7/ft), (18) 

where ^12)7=0 °c u {Aq/Tl) 1 ^ 2 eT Aa ^\ according to the general result de- 
scribed above. One can see that for half- integer 7/ft the interference of in- 
stanton and antiinstanton paths is destructive, so that at7 = ±|,±^,... 
the tunneling between two wells is completely suppressed. This effect is es- 
sentially topological because the topological charge appears in the answer: 
the contribution of configurations with different topological charge is dif- 
ferent. The same result can be obtained directly by solving the Schrodinger 
equation with the Hamiltonian (16): for half- integer 7/ft it can be mapped 
to the Mathieu equation with antiperiodic boundary conditions, and the 



corresponding energy levels are known to be doubly degenerate [25], which 
also means absence of tunneling. 

3. Field-Theoretical Description of a Small Ferromagnetic Par- 
ticle 

In this section we consider the basic technique of field-theoretical descrip- 
tion for spin systems on the simplest example, namely a nanoparticle of a 
ferromagnetic material. Such an object may be viewed as a zero-dimensional 
magnetic system, because at very low temperature all spins in the particle 
can be considered as pointing in the same direction. 

It is worthwhile to consider first the dynamics of a single spin S. In 
order to obtain the effective Lagrangian describing the spin dynamics, it 
is convenient to use a coherent state path-integral approach (see, e.g., the 
excellent textbook by Fradkin [26]). Let us introduce a set of generalized 
coherent states [27] 

|n) = exp{i0(n x l)S}\m = S) (19) 

parameterized by the unit vector n(0,(p). Here z is a unit vector pointing 
along the z axis, and \m) denotes a spin-S 1 state with S z = m. They form a 
non-orthogonal 'overcomplete' basis so that the following property, usually 
called a resolution of unity, holds: 

Jvn\n)(n\ = l, (20) 

another useful property is that quantum average of S on those coherent 
states is the same as of classical vector of length S: 

(n\S\n) = Sn . 

In case of S=l/2 those coherent states have a very simple form and are 
general single-spin wavef unctions: 

|n) = cos(0/2)| T) + sin(0/2)e i *'| |) . 

We again start from the formula for propagator (1) which is essentially 
a definition of the effective Lagrangian. Slicing the time interval [0; to] into 
infinitely small pieces At = to/N, and successively using the identity (20), 
one can rewrite this propagator in n-representation as 

Pab = lim / dn^dni - ■ ■ dn^{A\nQ){n^\B) 

N—*oo J 

x U(n k \e- iSAt/h \fi k+1 ) . (21) 

k=0 



Passing to the function n(t) of the continuum variable t, one ends up 
with the coherent state path integral (1) where the action A is determined 
by the effective Lagrangian 

L c s = -ih{(d t n\n) - (n\d t n)} - (n\H\n) . (22) 

It can be shown that the dynamical part of this Lagrangian has the 
form 

hS(l -cos 6>)^| ; (23) 

for arbitrary S this calculation requires some algebra, but for the simplest 
case S = \ it is straightforward. The expression (23) is nothing but the 
Berry phase [28] for adiabatic motion of a single spin. 

It should be remarked that the presence of the full derivative term 
hS(dip/dt) is rather nontrivial and allows one to capture subtle differences 
between integer and half-integer spins, as we will see below. For example, 
consider a single spin S in some crystal-field potential, with the effective 
Hamiltonian 

H = KS 2 Z - K'S 2 X , (24) 

where K, K' > and the easy-plane anisotropy K is much stronger than 
the in-plane anisotropy K'. The Lagrangian is 

L = HS(1 - cos 9) ^ - KS 2 cos 2 9 - K'S 2 sin 2 9 cos 2 <p . (25) 

There are two equivalent classical minima of the potential at 9 = ^, p = 
and 9 = ^, <p = it. Paths with 9 ~ ir/2 make the main contribution into 
the tunneling amplitude, so that we can approximately set 9 = | + 

C 1, and expand in i9 up to quadratic terms in the Lagrangian; in the 
term proportional to i9 2 the K' contribution may be neglected as small 
comparing to the contribution of K. After that, the "slave" variable f? can 
be excluded from the Lagrangian ("integrated out" of the path integral) 
because the corresponding equation of motion SL/Si} = allows to express 
i? through (p explicitly: 

Substituting this solution into the original Lagrangian (25), one obtains the 
effective Lagrangian depending on (p only: 



We see that we end up with the Lagrangian of a particle on a circle from 
the previous section, with the topological term 7 = hS. For each path 



where <p changes from to ir there is a corresponding antiinstanton path 
with ip changing from to — tt, and those paths contribute to the tunnel- 
ing amplitude with phase factors e lwS and e~ tnS . For half-integer S those 
contributions precisely cancel each other, making the tunneling impossi- 
ble. This is exactly in line with the well-known Kramers theorem, which 
states that in absence of external magnetic field all energy levels of a sys- 
tem with half-integer total spin should be twofold degenerate. One can also 
straightforwardly check that for a single spin in magnetic field, i.e., for 
H = gusHSz, the correct energy levels can be obtained only with the full 
derivative term taken into account. 

Now we are prepared enough, finally, to consider the problem of tun- 
neling in a small ferromagnetic particle consisting of N spin-S 1 spins. If we 
assume that ferromagnetic exchange interaction is so strong that we may 
consider all spins as having the same direction, then we come to the "giant 
spin" model where the entire particle is described as a quantum-mechanical 
("zero-dimensional") system with only two degrees of freedom 9 and p. In 
fact, we should postulate that in our path integral, when integrating over 
the coherent state configurations n^Li the main contribution comes 

from the subspace with all iV vectors Hi replaced by the same vector n(9, ip), 
and we take into account only configurations from this subspace. Assuming 
that the crystal- field anisotropy has the form (24), we come to essentially 
the same effective Lagrangian (27), and the only difference is that Eq. (27) 
should now be multiplied by the total number of spins N. The tunnel split- 
ting of the ground state level, according to Eq. (18), is given by 

AE = C{NS z ) 1,2 {KK' z ) l ' A \ cos{kNS)\ exp [-N S{2K' / K) 1 ' 2 ] , (28) 

where C is a numerical constant of the order of 1. A remarkable property of 
the result (28) is that presence of a large number ./V in the exponent can be 
to some extent compensated by smallness of the ratio K' /K. However, when 
the in-plane anisotropy K' — > 0, the splitting vanishes (this reflects the fact 
that in uniaxial case tunneling is impossible because of the conservation of 
the corresponding projection of the total spin; the same is true for K — > 0). 
Another remarkable feature is that for half-integer S the finite splitting 
can be observed only in particles with even number of spins N; since in any 
statistical ensemble N fluctuates a bit, this roughly means that only one 
half of all particles gives nonzero contribution. 

Statistical fluctuations of N have another, more painful consequence: 
since ./V stays in the exponent, even small fluctuations of the total number 
of spins in the particle lead to large fluctuations of the splitting. Moreover, 
since ./V scales as the third power of the linear size L, small fluctuations of L 
will be considerably enhanced in N. This may be crucial if one tries to detect 
the splitting by means of some resonance technique: the initially weak signal 



would be even more weakened by the strong broadening of the resonance 
peak. Actually, many factors can prevent one from observing the tunneling 
resonance, e.g., relaxation, temperature effects, etc. Here we will not at all 
touch the problem of relaxation because of its complexity; instead of that 
we refer the interested reader to the review [29]. Taking into account the 
finite temperature effects is also nontrivial, particularly because it requires 
changing the procedure of taking averages in the path integral: statistical 
averages should be taken simultaneously with quantum-mechanical ones. 
Roughly (and without taking into account the temperature dependence of 
relaxation mechanisms) the effects of finite temperature can be estimated 
with the help of the concept of a characteristic temperature T c below which 
the effects of quantum tunneling prevail over thermal transitions. Rough 
estimate for T c is obtained from the comparison of the relative strength of 
two exponential factors: thermal exponent e~^ u ^ T and tunneling exponent 
e ~Ao/h^ w ] iere £jj j s the- height of barrier separating two equivalent states 
and Aq is the corresponding instanton action, then T c = (HAU/Aq). It is 
easy to see that for the ferromagnetic particle problem considered above 

T FM = SiKK'12) 1 ' 2 , (29) 

i.e. the temperature of crossover from classical to quantum transitions is in 
this case rather small since it is determined by weak (relativistic) anisotropy 
interaction constants; for typical anisotropy values Tfm is about 0.1 K. 

4. Quantum Tunneling in a Small Antiferromagnetic Particle 1 

4.1. CONTINUUM FIELD MODEL OF ANTIFERROMAGNET 

The problem of continuum field description of antiferromagnet (AFM) is 
more complicated but also much more interesting than a similar prob- 
lem for ferromagnet. Antiferromagnet contains at least two different "sub- 
lattices" whose magnetizations compensate each other in the equilibrium 
state. Thus, when choosing the coherent state wavefunction in the form 

= Yii \fii) as described above, one cannot any more consider 
"smooth" function of the lattice site i. Let us adopt the simplest two- 
sublattice model which, despite the fact that it may be inadequate for 
a specific material, still allows one to demonstrate the essential physics 
of antiferromagnetism. We assume that there are two equivalent sublat- 
tices with magnetizations M\(f) and M2(f), \M\\ = \Mi = Mq. Then, 
when passing to the continuum limit, one has to introduce smooth fields 
fh = (Mi + M 2 )/2M and F= (Mi - M 2 )/2M describing net magnetiza- 
tion and sublattice magnetization, respectively. They satisfy the constraints 



Subsection 4.2.2 was written together with Vadim Kireev. 



ml = 0, rh 2 + r = 1, and we further assume that |m| <C |/|. The energy of 
AFM W = (H) then can be expressed as a functional of rh and I: 

W[rh,l]=M 2 j dV ^5m 2 + ±a(Vr) 2 + w a (l)-^(m-H)} . (30) 

Here the phenomenological constants 5 and a describe homogeneous and 
inhomogeneous exchange, respectively, H is the external magnetic field, 
g is the Lande factor, the function w a describes the energy of magnetic 
anisotropy, and we use the notation (V/) 2 = J2i{dl / dxi) 2 . The magnitude of 
sublattice magnetization Mq = grigS/vo, where \xb is the Bohr magneton, 
S is the spin of a magnetic ion, and vq is the volume of the magnetic 
elementary cell. 

As we learned from the previous section, the correct Lagrangian, suit- 
able for the quantum-mechanical treatment, has the form 

L = £>S{(1 - cos 6 U )^ + (1 - cos^)^} - W[mJ] , (31) 

i ' 

where the angular variables (6u, tpu) and (#2«, ^>2i) determine the unit vec- 
tors describing the orientation of spins in first and second sublattice, respec- 
tively. Note that we have kept intact the summation sign in the dynamical 
part of (31): the reason is that the explicit expression for the Berry phase 
in the continuum limit strongly depends on the details of the magnetic ele- 
mentary cell structure (which dictates the correct definition of fh and I and 
the procedure of passing to the continuum limit). 

Under the assumption that \ fh\ <C |Z|, the magnetization rh can be ex- 
cluded from the Lagrangian (31), and one obtains the effective Lagrangian 
depending only on I; after that step I can be regarded as a unit vector, 
P = 1. 

For example, in the simplest case of an antiferromagnet with only two 
(equivalent) atoms in elementary magnetic cell the dynamic part of the 
Lagrangian (31) can be written as 

J dV2hSm ■ (fx dl/dt) , (32) 

and the density of the effective Lagrangian takes the form 

^{^gy->i 2 -^+H r # <33) 

where w a is the anisotropy energy renormalized by the magnetic field, 

w a = w a + J ^(t-H) 2 , (34) 



7 = gus/fi is the gyromagnetic ratio, and c = ^jMo(aS) 1 ^ 2 is the limiting 
velocity of spin waves. Using general phenomenological arguments, one can 
show [30] that in case of arbitrary collinear antiferromagnet the Lagrangian 
should have the form similar to (33). 

Other, more complicated interactions can be present in Eq. (30). In 
some AFM materials (which are, strictly speaking, weak ferromagnets) the 
so-called Dzyaloshinskii-Moriya (DM) interaction is possible. It can be de- 
scribed by including the term D^rail^ under the integration sign in into 
(30), where is some tensor (which is not necessarily symmetric or anti- 
symmetric) . The origin of the DM interaction is rather nontrivial, and there 
is a number of "selection rules" excluding the possibility of its existence, 
particularly the DM interaction cannot exist (i) if there is an inversion cen- 
ter interchanging sublattices; (ii) if there is a translation which interchanges 
sublattices, i.e. if the magnetic elementary cell is larger than the elementary 
cell of the original crystal lattice. It can be shown [31] that presence of the 
DM interaction can be taken into account by the substitution 

H i-> H = H - ^M D (35) 

in the Lagrangian (33), where the components of vector D are defined as 
A = D ik l k . 

If there exists a sublattice-interchanging inversion center, another in- 
variant may be present in (30), namely Hi(m ■ dl/dxi) (here /ij are certain 
exchange constants). It is very important for the physics of AFM in one 
dimension, as we will see later. 



4.2. SPIN TUNNELING IN ANTIFERRO MAGNETIC NANOPARTICLE 

In case of a small particle one can consider m and I as being uniform 
throughout the particle, i.e. as not having any space dependence. Then, 
the Lagrangian (33) takes the form 

L = ^^\6 2 + sin 2 6ip 2 + 2^6 (H y cos H x sin (36) 

+ 2799 [H z sin 2 9 — sin 9 cos 6(H y sin ip + H x cos tp)] j — MQW a , 

where ./V is the total number of magnetic elementary cells in the particle, 
H e = 5Mq/2 is the exchange field, the dot denotes differentiation with 
respect to time, and we used angular variables for the vector /, 

l z = cos 9, l x + il y = sin 9 e 19 . 

There is another possible effect, typical only for antiferromagnetic particles: 
due to the boundary (surface) effects, the number of spins in two sublattices 



can differ from each other. In that case the Lagrangian (36) will contain 
the additional term 

hvS(l-coB0)(p, (37) 

which is essentially the Berry phase of v non-compensated spins. Such a 
sublattice decompensation in fact should be present in any ensemble of 
nanoparticles, so that v has certain statistical variation. 

The full Lagrangian (36) is rather complicated, and for the sake of clarity 
we will consider separately the effects of field and DM interaction. 

4.2.1. Tunneling in presence of external magnetic field 
Consider a small AFM particle with easy-axis anisotropy 

Wa = \f3(l 2 y + l 2 z ) 

in external magnetic field H perpendicular to the easy axis. Then the Eu- 
clidean action takes the form 

kNS f , f fd9\ 2 . Onf dl P\ 2 n - rr ■ 2 n^P 

- -iw; +sm "U) 

+ uj 2 [sin 2 9 sin 2 ip + (1 + j 2 H 2 /u^) cos 2 e\ j (38) 

dip 



+ 



ihvS f dr{\- cos 9)^- 
J dr 



where r is the imaginary time, and ojq = ^jMo(5P) 1 ^ 2 is the characteristic 
magnon frequency {%ojq is the magnon gap). 

There are two equivalent states A and B with opposite direction of I 
along the easy axis Ox, and obviously the most preferable instanton path 
is given by 9 = 7r/2, ip = <p(r). The instanton solution for (p is the same as 
in case of particle on a circle, and one-instanton action is 

4^ = ( 4cj o ± 2iri~fH) ± iirvS , (39) 

where ± signs correspond to instantons and antiinstantons. Thus, the tun- 
neling amplitude Pab is proportional to 

(^£) 1/2 exp (_^} cos{m , S + 2 ^ S(ff/ft)) , (40) 



and the corresponding magnitude of tunneling level splitting (proportional 
to | Pas |) oscillates with the period AH = (H e /2NS) when changing the 



external field. This period AH may be rather small, for typical values of 
the exchange field H e ~ 10 6 Oe and the number of spins in the particle 
N ~ 10 3 -7- 10 4 one obtains AH ~ 10 2 -=- 10 3 Oe. The effects of this type 
were studied in [32, 33]. 

The result (40) illustrates also another remarkable feature: in any ex- 
periment probing the response of the ensemble of AFM nanoparticles at 
each H there must be only one possible value of splitting (i.e., only one 
peak in the low-frequency response) when the spin of magnetic ions S is 
integer; but if S is half-integer then, since in any ensemble v arbitrarily 
takes even and odd values, for approximately one half of all particles the 
phase of cosine in (40) is shifted by 7r/2, and there should be two peaks at 
each H. 

It is worthwhile to note that the real part of the one-instanton action, 
which enters the exponent in (40), is proportional to (K/J) 1 / 2 (where J 
and K are the exchange and anisotropy constants) while the corresponding 
quantity for ferromagnet, according to (28), does not contain the exchange 
constant and is determined by the rhombicity (K' / K) 1 / 2 . One may con- 
clude that tunneling in AFM particles is more easy than in FM; indeed, the 
characteristic crossover temperature below which quantum effects dominate 
over thermal ones, for antiferromagnets is 

T AFM oc S(KJ) 1 / 2 , (41) 

which is much greater than for ferromagnets [cf. Eq.(29)]; typically Tafm 
is about 1 v3 K. 



4.2.2. Tunneling in presence of the DM interaction 

Consider the same small AFM particle from the previous subsection, but 
imagine that the DM interaction in its simplest form is present, with the 
energy given by 

w d = d(m y l z - m z ly) . (42) 

Then the DM interaction leads to the contribution into the Lagrangian (36) 
of the form 

AL d = • 2H D — (sin 6 cos if) , (43) 

7 .He 

where Hr> = dM is the so-called Dzyaloshinskii field. This term will con- 
tribute to the imaginary part of the Euclidean action (38), and as a result 
the cosine in (40) will be modified as 

cos {ttvS + 2TiNS(H/H e ) + ANS{H D /H e )} . (44) 



Thus, presence of the DM interaction alone also leads to effective change 
of the Berry phase and lifts the degeneracy for odd v and half-integer S. 



5. Spin Tunneling in Topological Magnetic Nanostructures 2 



As we mentioned before, one of the most difficult experimental tasks when 
trying to detect the resonance on tunnel-sphtted levels in small particles is 
to prepare the ensemble of particles with very sharp size distribution: even 
small fluctuations of size lead to large fluctuations of the tunneling proba- 
bility since they contribute to the power of exponent. Preparing such an en- 
semble requires high technologies and involves considerable difficulties. One 
may think about some other, "natural" type of magnetic nanostructures to 
observe spin tunneling phenomena in. One nice solution, which have actu- 
ally been used in experiment, is to use biologically produced nanoparticles 
[14]- 

Another possible way, proposed in [17, 18, 20], is to use topological^ 
nontrivial magnetic structures: kinks in quasi-lD materials, vortices and 
disclinations in 2D, etc. Such objects have required mesoscopic scale (e.g., 
the thickness of a domain wall is usually about 100 lattice constants) and, 
since their shape is determined by the material constants, they are identical 
to a high extent (up to a possible inhomogeneity of the sample) . 

Here we consider several possible scenarios of tunneling in topological 
nanostructures and show that their use has a number of advantages. 

5.1. TUNNELING IN A KINK OF ID ANTIFERROMAGNET 

Consider a one-dimensional two-sublattice antiferromagnet with rhombic 
anisotropy described by the Hamiltonian 

H = JJ2SA +1 + £ [K^Sff + K 2 {Sff] , (45) 

i i 

where i labels sites of the spin chain with the lattice constant a, K\ > K 2 > 
are the anisotropy constants (so that Oz is the difficult axis and Ox is the 
easy axis), and J is the exchange constant. For passing to the continuum 
field description one may introduce vectors m and I as rhk = (n2fc+i+^2fc)/2 
and Ik = (n2k+i — ^2fc)/2, where n are the unit vectors describing the 
direction of spins (the parameters of the corresponding coherent states, see 
the discussion in Sect. 4.1) above. These fields live on the lattice with the 
double spacing 2a, and it is easy to see that the energy functional W = (H) 
contains the term fn ■ d x l. Using the equation SL/Sm = 0, one may express 
fh through I and its derivatives and exclude it from the Lagrangian. The 
effective Lagrangian takes the following form: 

L cS = J g |^ (d t lf - JS 2 a\dJf - K.SHl - K 2 S 2 l^ + L top , (46) 

2 Subsection 5.2 was written together with Vadim Kireev. 



which represents a (l+l)-dimensional nonlinear a-model with the so-called 
topological term 

Lto P = \hS j ' dx f ■ (djx d t l) . (47) 

It is easy to trace the origin of this term: because of presence oifn-d x l in the 
energy, the expression for rh contains d x l which after the substitution into 
the Berry phase (23) yields the topological term. In agreement with gen- 
eral phenomenological result (33), the Lagrangian (46) is Lorentz- invariant, 
with the limiting velocity c = 2JSa/h. 

A stable kink solution corresponds to rotation of vector I in the easy 
plane (xy): 

l x =a'tanh(x/A), l y = J l z = , (48) 

cosh(x/A) 

where A = a(J j K^ 1 ! 2 is the characteristic kink thickness, and the quan- 
tities a and a' may take the values ±1. The topological charge of the kink 
a 1 is determined by the boundary conditions and cannot change in any 
thermal or tunneling processes. The situation is different with the quantity 
a which determines the sign of I projection onto the "intermediate" axis 
Oy. Two states with a = ±1 are energetically equivalent; change of a is 
not forbidden by any conservation laws and describes the reorientation of 
the macroscopic number of spins N ~ A/a » 1 "inside" a kink, typically 
N ~ 70 -=- 100. 

Again, tunneling between the kink states with a = ±1 can be studied 
using the instanton formalism. In contrast to the case of a nanoparticle, 
here the tunneling between two inhomogeneous states takes place, so that 
nontrivial space-time instantons come into play. The instanton solution 
Io(x,t) is now two-dimensional and has the following properties (see Fig. 
2): 

l x — > ±cr' at x — > ±oo 



l y — ► =Fc at x = 0, r — > ±oo (49) 
lz = V = ±1 at x = 0, r = . 

Along any closed contour around the instanton center in the Euclidean 
plane vector I rotates through the angle 2irv in the easy plane (xy), where 
v = aa' = ±1. Thus, the instanton configuration has the properties of a 
magnetic vortex and is characterized by two topological charges [34, 35]: 
vorticity v and polarization p. The instanton solution satisfies the equations 

V 2 ^ + sin0cos0[(l + psin 2 v3)/A 2 - (V(^) 2 ] =0, 

V • (sin 2 eV(p) - (p/A 2 ) sin 2 9 sin (p cos y = 0, (50) 




Figure 2. The structure of instanton solution for the problem of tunneling in a kink of a 
ID antiferromagnet. Arrows and circles denote projections of vector I on the easy plane 
(xy) and on the difficult axis Oz, respectively. Vector I forms the angle of about 45° with 
the easy axis Ox on thin solid curves, and with the difficult axis Oz on the circle (the 
circle radius is approximately ro). 



where we have introduced the angular variables l y +il z = sin 9e lip , l x = cos 9, 
p = [K\ — K2)/K2 is the rhombicity parameter, and V = (8/8x1,8/8x2) 
is the Euclidean gradient, (21,22) = (x,ct). 

Several important properties of the instanton can be obtained without 
using the explicit form of the solution. First of all, note that this instanton 
has two zero modes which correspond to shifting the position of its centre 
along the direction of r and 2 axes, respectively. The physical meaning 
of the first mode is the same as for ID instanton, and the second mode 
corresponds to moving the kink center in real space (the kink position 
in infinite ID magnet is not fixed in our continuum model); however, if 
the kink center is fixed due to some effects (e.g., because of pinning on 
the lattice, or by boundary conditions), so that the eigenfrequency of its 
oscillations is comparable with the characteristic magnon frequency, then 
only one zero-frequency mode is present. 

The Euclidean action Ae nc \ can be represented in the form 



Amci = -SHF + i2irShQ, where 



F=^J fx [(W) 2 + sin 2 8(V(p) 2 + cos 2 9 (1 + psin 2 <p)] 

Q = — J (fx £ a p sin OdaOdpip . (51) 

Imaginary part of the Euclidean action is in this case completely determined 
by the topological term L top . The word "topological" becomes now clear, 
because Q is the homotopical index of mapping of the (x\,X2) plane onto 
the sphere P = 1 (the Pontryagin index, or the winding number). For 
uniform boundary conditions at infinity in the {x\,X2) plane Q can take 
only integer values, but in our case Q = —pv/2 = ±^ is half- integer, which 
is typical for vortices (see, e.g., [34, 35]). For a kink with given a' there are 
two instanton solutions with the same vorticity v and different polarizations 
p. Thus, the tunneling amplitude is proportional to cos(ttS') and vanishes 
when the spin S of magnetic ions is half-integer. However, the degeneracy 
can be lifted in presence of external magnetic field or the DM interaction, 
as we will see below. 

We are not able to construct the exact solution of Eqs. (50), but the es- 
timate of the tunneling amplitude in various limiting cases can be obtained 
from approximate arguments. For p <C 1 the characteristic space scale of ip 
variation A/^/p is much greater than the kink thickness A, and the prob- 
lem can be mapped to one with a finite number of degrees of freedom (one 
may introduce the variable <ft having the meaning of the angle of deviation 
out of the easy plane "inside a kink" , so that the instanton solution can be 
seeked in the form cf) = 4>(t)), then it is easy to obtain [36] 

F ~ V /2 at p < 1 . (52) 

In the opposite limiting case /)> 1 one again has two different length 
scales: the kink thickness A and the "core" radius tq = A(K2/Ki) 1 / 2 , ro <C 
A. For r <C A all interactions except the exchange one can be neglected, 
and one may use the "isotropic" vortex solution 

e = e (r), ip = vx, v = ±i, 

d 2 9o ( 1 ^ 2 \ 

- d -^+ ( A2 - 72 J +sin# cos# = 0, (53) 

where r = (x 2 + x 2 .) 1 / 2 , x = arctan^/xi) are polar coordinates in the 
(xi, X2) plane. For r > ro, i.e., far outside the core, one can approximately 
assume that 

0=\, VV= ^sh^. (54) 

Within a wide range of r (for r <C r <C A) the solutions (53) and (54) can 
be regarded as coinciding, and the integrand in F is proportional to 1/r 2 . 



Then, one may divide the integration domain into two parts: r < R 
and r > R, where R is arbitrary in between ro and A. For r < R the 
solution (53) may be used, yielding F r< R = 7rln(£i?/ro) with ( ~ 4.2 [37]. 
For r > R, one can use a simple trial function approximately satisfying 
(54), e.g., 

x 2 1 xi 1 

CO sv3 = — sinp = , 55 

r cosh(x/A) r cosh(x/A) 

which yields -F r >R = 7rm(£'A/i?) with £ ~ 0.1. Summing up the two con- 
tributions, we obtain 

F~vrln(0.42A/r ) at p » 1 . (56) 
The tunnel splitting of the "ground state" level of the kink 

r oc hui (FS/2) n / 2 e- (FS/2) |$| , (57) 

where u)\ = 2S(JK 2 p) 1 / 2 is the frequency of the out-of-plane magnon local- 
ized at the kink, $ is the factor determined by the imaginary part of the 
Euclidean action [in the simplest model $ = cos(7r5)], and n is the number 
of zero modes which can be equal to 1 or 2 depending on whether the kink 
position is fixed, see above. It is easy to estimate the crossover temperature 
for the problem of tunneling in a kink, comparing the exponent in (57) with 
e~ u °/ T , where U ~ 2S 2 (^JK 1 - \J JK 2 ) is the barrier height; for p > 1 
(i.e., K\ » K-i) and n = 1 one obtains 

sjjK.y/ 2 

ln(K 1 /K 2 ) 

which is only logarithmically smaller than the corresponding temperature 
for a particle (41). 

Let us discuss now the behavior of the imaginary part of the Euclidean 
action in case of deviations from the simplest model (45) for which the 
tunneling is prohibited for half-integer S. The most simple observation 
is that in a spin chain with alternated exchange interaction, when along 
the chain the strength of exchange constant alternates as J\ J 2 J\J 2 ■ ■ ■, the 
topological term (47) acquires additional factor J\j J 2 (see, e.g., [38, 39]), 
which leads to $ = cos(ttSJi/J 2 ) and allows tunneling for half-integer S. 
Another way to lift the degeneracy at half-integer S is to "switch on" the 
DM interaction or external magnetic field. 

Consider the same model (45) with the addition of a magnetic field H 
applied in the easy plane (xy) . Presence of the field leads to the additional 



contribution to the imaginary part of A euc \ 



^CUCl 1 r ^CUCl ifoQ ; 

where n = H/H. The mixed product in (59) can be rewritten in angular 
variables as 

• „ /w n d<f . . 86 

— sin 8 cos (n x cos ip + n v sin </?) — h (n v cos ip — n x sin </?) — — . 

ax2 ' 0x2 

One may note that sin (9 and significantly differ from zero only in the 
vortex core, and thus the isotropic vortex solution (53) may be used for the 
calculation of Q' . After integration we obtain 

H A 

Q' = 2S— p(An x + vBn y ) , 

ti e a 

A= (dr I A) sin 0„ cos O , B= (dr/A)r(d9 /dr), (60) 
JO Jo 

where p and v, as earlier, denote the polarization and vorticity of the in- 
stanton solution, and A, B are numerical constants (recall that, according 
to (53), the isotropic solution #0 may depend only on r/A). After perform- 
ing the summation in p, u, and with the account taken of the contribution 
Q coming from the topological term, the factor <E> in (57) will be modified 
as 



$ ^ $ H = cos (2ASn x ——\ cos (irS + 2BSn v — — ) , (61) 
V H e aJ V H e aJ 

which means that for the given geometry only the field component per- 
pendicular to the easy axis lifts the degeneracy existing for half-integer S. 
Similarly to the case of a small AFM particle, the tunneling amplitude is an 
oscillating function of the external magnetic field H , but here the situation 
is more complicated because the period of oscillations depends on the field 
orientation. 



5.2. TUNNELING IN ANTIFERROMAGNETIC RINGS WITH ODD 
NUMBER OF SPINS 

Another example of a magnetic nanostructure is a ring formed by magnetic 
atoms; such rings may occur in a dislocation core of a 2D crystal as shown 
in Fig. 3, and the characteristic feature of this object is that the number of 
atoms in the ring is odd. Here we consider only antiferromagnetic rings. In 
terms of the vector I such a ring is a spin disclination. Let us assume that 



Figure 3. A "ring" around the core of dislocation in two-dimensional antiferromagnet. 
The dislocation is shown with a dashed line. 



the magnetic anisotropy is of the easy-plane type, and all spins lie in the 
(xy) plane, 

Si = (-I) 1 (e x cos Lpi + eysanpi) , 

where e x y are the unit vectors along x, y. Then there are two energetically 
equivalent states of the ring, with (fi = Xi/2 and ipi = — Xi/2, where Xi 1S 
the azimuthal coordinate of the i-th spin (let us assume that the ring is a 
circle of radius R). It is possible to construct the instanton solution which 
links the two states; in terms of I it can be written as 

l x = cos ^ , l y = sin ^ cos ip, l z = sin ^ sin ip , 

cosV' = ±tanh(oJoT), — -jlMo(j35) l l 2 . 

Calculation shows [22] that the tunneling amplitude is proportional to 

cos(^S) exp{-7rS12/A}, A = (a//3) 1/2 , (62) 

i.e., the probability of tunneling is sufficiently large if the radius of the ring 
is smaller than the characteristic thickness of the domain wall A (usually 
A ~ lOOA). Again, the tunneling is suppressed for half- integer S, and this 
can be changed with the help of external magnetic field. More detailed 
analysis [22] shows that the field H should be applied in the easy plane in 
order to lift the degeneracy, then the cosine in (62) will change into 

cos ( irS + Tr 2 S-^— — 

AH P a 



where a is of the order of the lattice constant. For weak fields the above 
expression describes just the Zeeman splitting of the ground state level of 
a ring (recall that due to the odd number of spins the ring always has an 
uncompensated total spin if S is half- integer) . 

5.3. TUNNELING IN A MAGNETIC VORTEX OF 2D 
ANTIFERRO MAGNET 

One more example of a magnetic topologically nontrivial structure is mag- 
netic vortex in quasi-2D easy-plane antiferromagnet. Consider the system 
described by the Hamiltonian 

H = jJ2Si-S j + Kj2(S!) 2 (63) 

where K > is the anisotropy constant, and Oz is the difficult axis. In 
terms of the angular variables for the antiferromagnetism vector, l z = cos 9, 
l x + il y = sin9e llp , a vortex corresponds to the solution 

= 0%{r), <p = v X + <Po, (64) 
0±(oo)=7r/2, 9+(0) = 0, 9 (0)=n, (65) 

where 9q satisfies the equation from the second line of Eq. (54), x+iy = re iX , 
and the solutions 9q and 9q have the same vorticity v but different polar- 
izations p = cos 9(0) = ±1. The vortex states with p = ±1 are energetically 
equivalent, and the transition between them corresponds to reorientation 
of a macroscopic number of spins N ~ (A/a) 2 , where A = a(K/AJ) 1 / 2 is 
the characteristic radius of the vortex core and a is the lattice constant. It 
is worthwhile to remark that such a transition would be forbidden in fer- 
romagnet because of the conservation of the z-projection of the total spin 
S z . 

The instanton solution l(x,y,r) linking two vortex configurations 9^ 
with v = 1 when the imaginary time r changes from — oo to +oo is schemat- 
ically shown in Fig. 4. In the 3D Euclidean space (x, y, r) it describes a 
topological configuration of the hedgehog type and has a singularity at the 
origin. Such a singularity means that in a small space region around the 
origin (roughly within the distance of about a) one has to take into account 
the change of magnitude of the sublattice magnetization: the length of the 
vector I has to change so that |/(0, 0,0) | = 0. In this case there are four 
zero modes, three of them correspond to translations along x, y, r, and 
the fourth one corresponds to changing the (fo angle. If the position and 
structure of the vortex are fixed by some additional interactions, only one 
zero mode is left. 




The Euclidean action derived from the Lagrangian of the a-model has 
the following form 



A E = JS 2 /dr/^x 1 1 ^ j 2 + (Vf) 2 - (V/™) 2 + - (Zf) 21 

(66) 

where ft ) describes the vortex solution (64) and c denotes the limiting 
velocity c = 2JSa/H. Away from the singularity (for p » a, p = (c 2 t 2 + 



x 2 + y 2 ) 1 / 2 ) the condition r = 1 holds, and the equations for 9, ip become 

V 2 # + sin#cos0[l/A 2 - (Vtp) 2 } = 0, 

V • (sin 2 OVcp) = . (67) 

In the region a <C p <C A this system has an exact centrally symmetric 
solution of the hedgehog type: 

ct y 
cos 9 = — , tan 09 = —. (68) 
p x 

It can be shown that the contribution of the singularity itself is small and 
can be neglected. Dividing the integration domain into two regions p < R 
and p > R, where i!« A, one can see that the contribution of the region 
of small distances p < R to the Euclidean action is given by 

A E [p <R}= Ait{JS 2 /c)R . (69) 

For estimating the contribution of the "large" distance region we use a 
variational procedure with the trial function of the form 

0(x, y, t) = vr/2 + F(cr) [vr/2 - 9 { +) (r)] , (70) 

where F(ct) is a "smeared step function:" F — > ±1 as r — > ±oo and the 
derivative of F is nonzero in the region of the thickness Ai around r = 0. A 
simple estimate shows that the resulting contribution of the region p > R 
is described by 

A E [p >R] = (2vr JS 2 /c) [Ci Ax ln(A/i?) + £ 2 Ax + 6A 2 /Ax] , (71) 

where £1,2,3 are numerical constants of the order of unity. Summing up (69) 
and (71) and minimizing Ae with respect to Ai and R, we find Ai ~ R ~ 
A. Thus, the total one-instanton Euclidean action may be estimated as 

A = 2tt£JS 2 A/c = ZTrhSA/a , (72) 

where £ ~ 1. 

Demanding that the tunneling exponent is not too large, e.g., Aq < 
20 -7- 30, we see that for S = 5/2 this means A/a < 3 4- 4, which is rather 
tight; the continuum field approach we used here formally requires A S> a, 
but in practice it is still applicable for A/a ~ 2 4-3 [40]. The crossover 
temperature T c ~ S(JK) 1 / 2 is not small since it is proportional to \fl . 

6. Summary, and What is left under the carpet. 

Let us mention briefly the problems which are closely related to the topic 
of this paper but were left out of discussion, and also those problems which 
are not clear at present, to our opinion. 



First of all we would like to remark that we did not touch at all micro- 
scopic essentially quantum effects in magnets, e.g., predicted by Haldane 
destruction of (quasi) long-range order in ID antiferromagnets with integer 
spin S caused by quantum fluctuations. Effects of quantum interference 
are also important for this phenomenon, and its existence is determined 
by the presence of topological term in the Lagrangian of antiferromagnet 
(see, e.g., the reviews [34, 39]). For small S and weak anisotropy the ground 
state of f D antiferromagnetic system can differ drastically from its classical 
prototype; e.g., the ground state of a S = 1 AFM ring is not sensitive to 
whether the number of spins is odd or even and is always unique, and the 
ground state of a S = \ AFM ring with odd number of spins is fourfold 
degenerate [22]. 

We also did not consider the contribution of tunneling-generated in- 
ternal soliton modes to the thermodynamics and response functions of ID 
antiferromagnets, which can lead to interesting effects (see [36, 35, 19, 41]). 

Another problem which was ignored in our consideration is the role of 
relaxation and thermal fluctuations of different origin. Even at low temper- 
ature the interaction of spins with other crystal subsystems (lattice, nuclear 
spins, etc.) may be very important, see [42, 43, 44]. It is clear that stochas- 
tic influence on the dynamics of magnetization from thermal fluctuations 
leads to decoherence and suppresses coherent tunneling. Description of this 
fundamental problem in any detail goes far beyond the scope of the present 
lecture, and we refer the reader to the review by Caldeira and Leggett [29] . 

One more problem which is unclear from our point of view is a justi- 
fication of considering all spins in a small particle as moving coherently 
(a "giant spin" approximation usually used in treating the MQT problems 
and also adopted in the present paper). In fact, the only justification of 
this approximation is energetical: if the particle size is much smaller than 
the characteristic domain wall thickness, any inhomogeneous perturbation 
costs much energy. On the other hand, for the Hamiltonian (24) neither 

S 2 nor S z are good quantum numbers, which means presence of magnons 
(deviations from collinear order) in the ground state. 

Our lecture was devoted first of all to the fundamental aspects of MQT 
considered as a beautiful physical phenomenon which is rather difficult to 
observe. But technological development can lead to the situation when this 
phenomenon will become practically important. The present tendency of 
increasing the density of recording in the development of information stor- 
age devices means decrease of the elementary magnetic scale corresponding 
to one bit of information, and one may expect that quantum effects will 
determine the "natural limit" of miniaturization in future. 
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